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Magnetic flares and induced oscillations of magnetars (supermagnetized neutron stars) are promis- 
ing sources of gravitational waves (GWs). We suggest that the GW emission, if any, would last longer 
than the observed x-ray quasiperiodic oscillations (X-QPOs), calling for the longer-term GW anal- 
yses lasting a day to months, compared to than current searches lasting. Like the pulsar timing, 
the oscillation frequency would also evolve with time because of the decay or reconfiguration of the 
magnetic field, which is crucial for the GW detection. With the observed GW frequency and its 
time-derivatives, we can probe the interior magnetic field strength of ~ 10^'' G and its evolution 
to open a new GW asteroseismology with the next generation interferometers like advanced LIGO, 
advanced VIRGO, LCGT and ET. 

PACS numbers: 95.85.Sz,97.60. Jd 
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In the upcoming years, the gravitational wave (GW) 
astronomy will be started by the 2nd generation GW de- 
tectors [l|-y|, such as advanced LIGO, advanced VIRGO 
and LCGT, and the 3rd generation ones like ET [ij in the 
10 Hz-kHz band. These GW interferometers will bring 
about a greater synergy among multimessenger (electro- 
magnetic, neutrino, cosmic ray and GW) signals. One 
of the most important targets is the neutron star (NS) 
oscillations. We can infer the NS interior from eigenfre- 
quencies of the waveform, and could open gravitational 
wave asteroseismology 5, 6]. 

Highly magnetized NS's with lO^'^'^G so-called 
magnetars [utM promising candidates for the GW 
asteroseismology. They are observed as soft gamma-ray 
repeaters (SGRs) and anomalous x-ray pulsars (AXPs) 
in our galaxy and the LMC, whose emission is thought to 
be powered by the dissipation of magnetic energy. SGRs 
emit repeated bursts and sometimes gian t flares (GFs) 
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source is the magnetic field, not 



the GW-scopy of the internal 
The magnetar GFs would excite the seismic oscillations 



with GWs |14| - |16| . The GFs are thought to be produced 
by a release of the accumulated magnetic stress in the 
crust or inner region. A part of the GF energy will tap 
the stellar oscillations, especially polar modes that vary 
the moment of inertia and emit GWs [l^ • At least the 
reduction of the magnetic stress would change the stellar 
deformation and begin polar oscillations around a new 
equilibrium. 

Actually, the quasiperiodic oscillations with ^ 100 
Hz have been discovered in the x-ray tails of GFs (X- 
QPOs) [3] and recently activated the searches of GWs 
from magnetar GFs. The GW energy within the X-QPO 
duration ~ 200 sec has been limited below ^ 10% of 
the electromagnetic energy X-QPOs can be inter- 

preted by the axial oscillations propagating to the mag- 
netospheric emission via magnetic fields 18l420| . 



In this Letter, we suggest that the GWs from magne- 
tar GFs would last much longer than the X-QPOs if the 
GW frequencies are close to the X-QPO frequencies [see 
Eq.®]. We show that the long-term GW analyses from 
a day to months are necessary to detect the GW, even 
if the GW energy is comparable to the electromagnetic 
energy. The current searches of < 200 sec are not long 
enough; this is because the GW luminosity is propor- 
tional to (frequency)^. Thus, by fixing the GW energy 
budgets, ^ lOOHz GWs could persist up to 10^ times 
longer than ^ kHz GWs of typical fluid modes like p- 
modes. 

In most previous studies, the fluid modes were em- 
ployed for GWs. However, ^ 100 Hz oscillations rather 
than ^ kHz fluid modes will be excited more effectively 
in the magnetar GFs because the initial spike of GF had 
a rise time '-^ 0.01 sec (lol-[l^ and its inverse ^ 100 Hz 
is most resonant. Actually X-QPOs have this frequency. 
This is also about the internal Alfven crossing frequency 
[see Eq.Q], which is one of the reasons to believe the 
magnetic origin of GFs [8!|. 

So, we propose that GWs from magnetar GFs, if 
any, are largely produced by the polar Alfven oscilla- 
tions, which last longer than the axial-type X-QPOs. 
This opens an interesting possibility to directly measure 
the internal magnetic field from the observed frequency 
faOcB since the restoring force is magnetic and the polar 
Alfven modes have discrete frequencies 17 1. 

The long-term GW emission also makes it possible to 
directly measure the magnetic field decay or reconfigu- 
ration inside the NS, because the mode frequency / de- 
pends on B that evolves with time. The Taylor expanded 
phase of the GW, 



*(t) = V^o + 27r <; fit - to) - ^/(i - to)' 



(1) 



has the first frequency derivative (IFD) term and also 
the higher frequency derivative terms (2FD and so on) 
like the pulsar timing search. The FD terms affect the 
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signal to noise ratio (S/N) even if we assume the mini- 
mum magnetic field decay required to supply energy for 
the observed emission from magnetars [see Eq.®]. 

In our picture, GWs are produced by the polar oscil- 
lations while X-QPOs are caused by the axial ones. This 
is implied by the amplitude of X-QPOs ^0.1 [18], which 
is much larger than the expected polar oscillation am- 
plitude < 10"* [see Eq. ©]. Thus, the GW searches 
beyond the X-QPO duration are meaningful. 

To clarify the above statement, we first make an 
order-of-magnitude estimate before the detailed analy- 
sis. We take the fiducial minimum magnetic field in- 
side the NS as B = lO^^G because the total magnetic 
energy ~ (B'^/Stt) ■ (AttR^/S) ~ 10"^^ erg should sup- 
ply - lO^^erg X lO^yr/lO^yr lO^^erg for the GF ac- 
tivities once per century (with three GFs out of ~ 10 
sources within 30 years) during the lifetime of magnetars, 
Tmag ~ 10'*yr, and the ionospheric observations suggest 
that ~ 9 times energy is released outside the X-7 ray 
band [211]. 

When the interior magnetic fields are deformed by the 
GFs, the induced difference in the moment of inertia of 
the star e is given by 



AI 
T 



^ [By St:) ■ [A^R" 1^) 



3 X 10-^Bi6^ 



(2) 

where / « 2MR^/5 is the moment of inertia of the 
magnetar, Biq = (i?/10^^G) is the interior magnetic 
field, and M, R are the mass and the radius, respec- 
tively [13, [11]. Hereafter, we fix i? = lO^cm and M = 



IAMq. The oscillation energy or the GW energy is about 
the gravitational energy shift caused by the stellar defor- 
mation, 



Egw ~ (eVs) • [GM^/R) < 8 x lO^^Big* 



[ergj 



(3) 



The oscillation frequencies are determined by the Alfven 
sound crossing, 

1 B 
R ' 



IOOB16 [Hz 



(4) 



where Va is the Alfven velocity, and p w M/ {AttR^ /3 
the mass density. If these oscillations are polar mode 
in which the moment of inertia changes, GWs are emit 
ted. The luminosity of the GWs are estimated to be 



ij ^ ij 



< 2 x W^B 
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(G/5c5) • (e/(27r/,)3)2 
erg/sec]. (5) 



Then the GW duration after the GF becomes 

TGw ~ Egw/Lgw 4 X lO^Bie"'^ [sec]. 



(6) 



The characteristic GW amplitude he [23[ is given by 
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TGW 



< 5 X 10-23B^g3/2^^^-l^ 



(7) 



where dig = d/lOkpc is the distance from the source. 
Therefore, if the GW energy is > lO'^'^erg (only a fraction 
of the GF energy ^ xq^''"^^ erg) and the data taking 
lasts for a day to months, the GWs from magnetar GFs 
within 10 kpc are detectable with enough S/N by the next 
generation detectors. The event rate will be ~ 1/lOyr, 
and more for weaker flares. We should note that the GW 
energy or the deformation e has a large uncertainty. The 
minimum deformation would be e ~ Cmax/lOO in Eq.Q 
since ~ 1% of magnetic energy is released in a GF. We 
can rescale the following results with S/N oc he oc e. 
Meanwhile, Eqs.Q and ([6]) do not depend on e, and thus 
are relatively robust. 

Since the magnetic fields of magnetars are considered 
to decay significantly in the lifetime T^ag ^ 10'*yr, the 
magnetic field are decaying at least at the rate 
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-B/t, 



mag 



-10^^16 [G/sec] 



This leads to the evolution of the Alfven frequency 



-IQ-^Bie [Hz/scc]. 



(8) 



(9) 



During the GW emission, the oscillation frequency 
changes by 



|A/a| ~ I/a X TGW\ 4 X lO-^Big-^Hz. 



(10) 



In principle, the frequency resolution of GW is given by 
the inverse of the observational time Tobs, 



Mob 



T, 



obs 



TGW 



"1 - 3 X lO^'^Bie^Hz (11) 



for long enough data taking. Therefore, we could resolve 
the GW frequency evolution to prove the magnetic field 
decay rate inside the magnetar for Biq < 3.0. In other 
words, the IFD term could be crucial for keeping S/N in 
the GW detection. If the magnetic field does not decay 
constantly in the GF phase, higher FD terms may be im- 
portant, where the long-term frequency-time analysis (2^ ] 
may be useful. 

Next, we shall move on to more details with the 
matched filtering Fisher analyses. We only consider the 
I = 2,m = mode as the magnetar oscillation because 
one-sided outflow has been confirmed by the radio obser- 
vation of the GF from SGR 1806-20 This suggests 
that a GF is a pointlike energy release from a NS, which 
would mainly trigger the I — l,m — plus I — 2,m — 
mode oscillations. We take m = since there is no spe- 
ciflc direction along the azimuthal angle and the rotation 
is slow (^ 0.1 Hz) compared with the Alfven oscillations. 
The Z = l,m = mode is the dipole oscillation which 
does not contribute to the GW emission. 

Taking a misalignment between the polar oscillation 
and the rotation axes as in FiglU the waveform of the 
observed plus (-I-) and cross (x) GWs are obtained from 
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the quadrupole formula, 
(I 



where 



(yl + i6)cos*(i)exp 



(12) 



+ sin^ e{i cos^ a - 1) + sin 261 sin 2a cos(27ri/t + 0) } (13) 



S = — ^ I — I {cos sin a sin 2(27ri'i + (^) 
V3 V / 

+ sin6'sin2Q!sin(27riyt + 0)} . (14) 

The stellar oscillation phase 4'(i) can be Taylor expanded 
as Eq.((T]). In this Letter, we stop the expansion at the 
IFD term, and add V'o, fa and fa as the parameters for 
the GW waveform. We neglect the proper motion of the 
magnetar. As a whole, the parameters are the strength 
of deformation e, the mean moment of inertia /, the dis- 
tance to the magnetar d, the rotational frequency of the 
magnetar u, the angle between the magnetar rotation and 
the oscillation axes a, the frequency of the oscillation /a, 
its evolution rate /a, the initial phase of the oscillation 
-00, the duration of the GW emission rciy, and the obser- 
vational polar and azimuthal angle 6 and 0. The position 
(0, (f) can be determined independently by the x-ray or 
radio observations. Here, /, d and e are completely de- 
generated. We fix / = lO'^^erg • sec^, d = lOkpc and 
leave e as a free parameter. According to Eqs.([2|),(|4l),(l6]) 
and ©, we set fa = lOOBigHz, = -lO-^BigHz/sec, 
TGW = 4 X 10*'Bi6~^sec and e = 3 x lO^^Sie^ (i.e., op- 
timal case). We also set a = 7r/4, v — O.lHz, 6 — 7r/2, 
<j) = -K /A and ip^ = 7r/4. 
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FIG. 1: The configuration of tlie magnetar and detectors. 



In the matched filtering analysis, we can compute the 
determination accuracy of these parameters using the 
Fisher matrix formalism with a waveform template in the 
frequency domain ft, (/) [26, 27|. The variance-covariance 



matrix of the parameter estimation error is given 
by the inverse of the Fisher information matrix Tij as 
(A7iA7j) — {T^^)ij. The Fisher matrix becomes 



4Re I 



a— +, X ' 



df dh*M)dhM) 

Snif) dj, dj, 



(15) 



where Snif) is the noise spectrum. The accuracies of the 
interested parameters are described as AB/B = A fa/ fa, 
AB/B = Afa/fa and AEgw/Eqw = 2Ae/e. The reso- 
lution for the position of the magnetar is 



An = 27rv/(A^2)(A02) - (A^A0)2, 
where /i = cos6'. The S/N is given by 

df 

I S 

a=+,x ■ 



(S/N)2=4 J2 I 



Snif) 
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(16) 



(17) 



We integrate the waveform in Eqs. ([T5)) and ([T7)) for the 
duration Tobs- We refer advanced LIGO (2nd gener- 
ation) and FT [4] (3rd generation) for a detector noise 
spectrum Snif)- 



TABLE I: S/N and parameter determination accuracy 
for the optimal case {e = tmax) 



detector 


2nd generation 


3rd generation 




1.0 


2.0 


1.0 


2.0 




Tobai^ tgw) 


4 months 


2 days 


4 months 


2 days 




S/N 


11 


23 


170 


350 




AB/B 


6.4 X 10~" 


1.0 X 10"** 


4.0 X 10"" 


6.5 X 10" 


10 


AB/B 


8.9 X 10~^ 


9.0 X 10"^ 


5.6 X lO"'^ 


5.8 X 10" 


-4 


Av/v 


4.2 X 10"** 


1.3 X lO"'^ 


2.7 X 10"® 


8.4 X 10" 


-8 


AEgw / Egw 


0.27 


0.14 


1.7 X 10"^ 


8.8 X 10" 


-3 


Atgw /tgw 


0.18 


0.090 


1.2 X 10"^ 


5.8 X 10" 


-3 


Aa/a 


0.030 


0.015 


1.9 X 10"^ 


9.8 X 10" 


-4 


An 


0.044 


0.010 


1.7 X lO""* 


4.2 X 10" 


-5 



TableU shows the S /N and the parameter determina- 
tion accuracy assuming enough data taking, Tobs ^ tgw- 
For the optimal case (e = Cmax), even 2nd generation de- 
tectors can detect the GWs with enough S/N. The de- 
tection of the GW is possible when e > O.lemax and e > 
0.01e„iax for 2nd and 3rd generation detectors, respec- 
tively. Thanks to the large Q value ~ /q x tqw ^ 10®, 
parameters related to the phase of the GW (i.e., fa, fa 
and v) could be determined accurately. When fixing B, 
errors scales as (S/N)"^ for ASl, and (S/N)"^ for other 
parameters. 

FigslHandOshow the dependence of the S/N and the 
accuracies on the data taking time Tobs- We can clearly 
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FIG. 2: The contour lines of S/N in the plane of the 
magnetic field of magnetars and the GW observation time. 
Solid (dashed) lines show S/N = 1, 10, 25, 45 (S/N = 10) using 
a 2nd (3rd) generation detector. B — 10^^ G is the minimum 
interior magnetic field required for magnetar activities. 
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FIG. 3: The parameter determination accuracies with respect 
to the GW data taking time. The solid, dotted and dotted- 
dash lines show \B\, Egw and AQ, respectively. Thin (thick) 
lines shows the results using 2nd (3rd) generation detectors. 
We set B = 1.5 X 10^*^0. 



find that the long-term data taking for a day (~ lO^sec) 
to months (~ lO^sec) is necessary for high S/N and pa- 
rameter accuracies, especially B. 

Since a NS has various oscillation modes, it is hard, at 
this stage, to answer which oscillation mode is mainly ex- 
cited at the GF. Also, we should consider the mode cou- 
pling between the polar Alfven mode and other modes. 
If the coupling is strong, the polar Alfven mode oscilla- 
tion cannot last so long as '-^ a day to months even if 



the GF energy are initially injected to the polar Alfven 
mode, and the dependence of Alfven mode frequencies 
on the magnetic field may become different from Eq.(|3|). 
The effects of the crust and superfluidity are also to be 
investigated in the future. 

The change in the rotational frequency v could be- 
come important as the IFD. Although the magnetar 
spindown rate is usually small ~ — 10~^^Hz/sec com- 
pared with Eq.([3]), the sudden increase in the GF phase 
Aiy/v ^ —10^^ has been reported in the August 27, 1998 
event from SGR 1900+14 |28|. Since v ^ fa, we could 
also measure i> independently. 

Finally, we should mention the strategic change in the 
analysis of the GW waveform. When discussing GWs 
from a single NS or black hole, the waveform usually 
has a short damping time < lOsec [2^ [s^, hence the 
evolution of the waveform due to the daily and yearly 
motion of the earth could be neglected. However in the 
case of the long-term oscillations with > a day, we have 
to take into account these effects as in the analysis of 
GW from a pulsar [3ll,|32|. 
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